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Abstract. We study the problem of determining whether a given frame 
is scalable, and when it is, understanding the set of all possible scalings. 
We show that for most frames this is a relatively simple task in that the 
frame is either not scalable or is scalable in a unique way, and to find 
this scaling we just have to solve a linear system. We also provide some 
insight into the set of all scalings when there is not a unique scaling. 
In particular, we show that this set is a convex polytope whose vertices 
correspond to minimal scalings. 



1. Introduction 

A collection of vectors {(fi}f =1 C C d is called a frame if there are positive 
numbers A < B < oo such that 



A\\x\\ 2 <^2\(x,^)\ 2 <B\ 



i=i 

for every x in C d . If we have A = B we say the frame is tight, and if 
A = B = 1 we say it is a Parseval frame. Given a frame {<£>i}" =1 we define 
the frame operator S : C d — > C d by 

n 

(i) Sx = y](x,<pi)<pj. 

It is easy to see that S is always positive, invertible, and Hermitian. Further- 
more, {<fii}2 = i is a Parseval frame if and only if S = 1^ (the identity operator 
on C d ). By a slight abuse of notation, given any set of vectors {y>i}E=i C C d 
we will refer to the operator defined in ([TJ as their frame operator, even if 
they do not form a frame (in this case S will not be invertible, but it will 
still be positive). If we have that \\<pi\\ = 1 for every i = 1, ...,n we say it is 
a unit norm frame. For more background on finite frames we refer to the 
book @|. 

A frame {<£>i}" =1 is said to be scalable if there exists a collection of scalars 
{vi}f =1 C C so that {vi(fi}™ =1 is a Parseval frame. In this case, we call the 
vector (|fi| 2 , \v n \ 2 ) G a scaling of {y>i}™ = x- Scalable frames have been 
studied previously in [5]. 
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We will work in the space Mdxd of all dx d Hermitian matrices. Note that 
this is a real vector space of dimension d 2 (it is not a space over the complex 
numbers since a Hermitian matrix multiplied by a complex scalar is no longer 
Hermitian). The inner product on this space is given by {S, T) = Trace(ST) 
and the norm induced by this inner product is the Froebenius norm, i.e., 
(S, S) = \\S\\ p. 

In what follows we will always consider frames in the complex space C d , 
however all of our results hold in the real space M. d as well. The only differ- 
ence is in this case we must replace the space Mdxd with its subspace §dxd 
consisting of all d x d real symmetric matrices, which is a real vector space 
of dimension d(d + l)/2. Thus, if one replaces H^xd with §>dxd and d 2 with 
d(d + l)/2 all of our results will hold for frames {(/?j}™ =1 C M. d and the same 
proofs will work. 

2. Scaling generic frames 
Consider the mapping from C d to Mdxd given by 

x h-> xx* . 

Note that xx* is the rank one projection onto spanjx} scaled by ||x|| 2 . xx* 
is called the outer product of x with itself. Also note that if x = \y for 
AgC then xx* = (Xy)(Xy)* = \X\ 2 yy*. 

Given a frame {^t}" =1 , in this setting we have that the frame operator is 
given by 

n 

s = ^2<Pi<p*, 

i=l 

so {y?i}™ =1 is scalable if and only if there exists a collection of positive 
scalars {u>j}™ =1 so that 

n 

^Wilfilfl = I d , 
i=l 

in this case {\A*Wi}iLi is a Parseval frame, and the vector (w\, 6 M" 

is the scaling. 

Before stating our first theorem we need one more definition. A subset 
Q C R n is called generic if there exists a polynomial p(x\, ■■■,x n ) such that 
Q c = {(xi, ...,x n ) S W 1 : p(x\, ...,x n ) =0}. It is a standard fact that generic 
sets are open, dense, and full measure. When we talk about a generic set in 
C d we mean that it is generic when we identify C d with M. 2d . 

Theorem 2.1. For a generic choice of vectors {^ij^i Q C d we have that 
span^jt/?*}^ = M d xd- 

Proof. First let {Ti}fLi be any basis for Mdxd- Since each T{ is Hermitian 
we can use the spectral theorem to get a decomposition T{ = Yll=i ^ij^ij 
where each is rank 1. So it follows that span{Pjj} = Mdxd and therefore 
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this set contains a basis of M d xd- Thus, we have constructed a basis of M^xd 
consisting only of rank 1 matrices. 

Now observe that for a given choice of vectors {(fi} d =i we have that 
sp&n{ifiip*} = Mjxd if and only if the determinant of the frame operator 
is nonzero (note that we are refering to the frame operator of {fi^pX^iLi as 
an operator on M^xd, not the frame operator of {<fi}f =1 as an operator on 
C ). But the determinant of the frame operator is a polynomial in the (real 
and imaginary parts) of the entries of the (fi's, and by the first paragraph 
we know that there is at least one choice for which this does not vanish, so 
we can conclude that for a generic choice it does not vanish. □ 

Corollary 2.2. If n < d 2 then for a generic choice of vectors {ifi}™ =1 C <C d 
we have that {^WiliLi is linearly independent. 

Given a frame {^j}™ =1 Q C d define the operator A : W 1 —> M-dxd by 

n 

Aw = y~]wnpi(p* 
i=i 

where w = (w\, w n ) T . To determine whether {(/?j}™ =1 is scalable boils 
down to finding a nonnegative solution to 

Aw = I d . 

In the generic case when {ifiip*}^ =1 is linearly independent, this system is 
guaranteed to have either no solution, or one unique solution. So if it either 
has a solution with a negative entry or has no solution we can conclude that 
this frame is not scalable, and if it has a nonnegative solution then it is 
scalable and this solution tells us the unique scalars to use. We summarize 
this in the following corollary: 

Corollary 2.3. Given frame {ipi}f =1 C C d such that {tpi(p*}f =1 is linearly 
independent in BL^xd, we can determine its scalability by solving the linear 
system 

(2) Aw = I d . 

Furthermore, in this case if it is scalable then it is scalable in a unique way. 

In particular, if n < d 2 then with probability 1, determining the scalability 
°/{ ( / 9 «}r=i i s equivalent to solving the linear system given in (|2|). 

3. Linearly dependent outer products 

In this section we will address the situation when {<Pi<p*}2=i is linearly 
dependent. The main problem here is that the system Aw = Id may have 
many solutions, and possibly none of them are nonnegative. In this section 
we will find it convenient to assume that \\<Pi\\ = 1 for every i = 1, n, note 
that we lose no generality by making this assumption. 
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Given a collection of vectors {xj}™ =1 C M d we define their affine span as 

n n 

aff{xj}" =1 := {^axi : ^ q = 1} 
i=l t=l 

and we say that {xj}" =1 is affinely independent if 

Xj aff {a^j-j^- 
for every j = 1, ...,n. We also define their convex hull as 

n n 

conv{xj}™ =1 := (^CiXi : a > 0, ^ Cj = 1}. 

i=l i=l 

We say a set "P C M. d is called a polytope if it is the convex hull of finitely 
many points. 

Proposition 3.1. Given a collection of unit norm vectors {9?j}™ =1 C C d 
we /iai>e that {ipi(f*}^ =1 is linearly independent if and only if it is affinely 
independent. 

Proof. Clearly linear independence always implies affine independence. So 
suppose that {¥W*}f = i is not linearly independent. Then we have an equa- 
tion of the form 

for some j. Also note that since \\(pi\\ = 1 it follows that (<pi<p*,Id) = 1 for 
every i = 1, ...,n. Therefore, we have 

1 = (<Pj <Pj i I d) = (^ WW** ' T d) 

Therefore {yw*}™ = i is not affinity independent. □ 

Proposition 3.2. A unit norm frame {</?j}" =1 C C d is scalable if and only 
\ld £ co^iVi^iliLi- Furthermore, if Xld G conv{<^j</?|}" =1 i/ien A = ^ 
and «/ XlILi w mP* = \h then J2i=i w i = 1- 
Proof. Suppose we have a scaling io so that 

i=l 

Then 

r( 

d = {Id, Id) = CyZ w ifi^iJd) 



i=l 

n n 

= Yl w i{ L Pi<-P*iJd) = Y W i- 
i=l i=l 
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Thus, Y17=l it = 1 an d since toj > for every i = l,...,n it follows that 
\h = Ya=i itPWi G conv { < ^j < ^*}r=i- The converse is obvious. 

The furthermore part follows from a similar argument. Suppose XI d = 
Yh=i Wi^fiLf* with Yh=i w i = L Trien 

n 

dX = (AJ d , Id) = y^Wj = 1. 

i=i 

Now suppose i/ d = X^iLi «W<Pi- Then 

n n 

1 = Wi(pi<p*,I d ) =^2wi. 

1=1 1=1 

□ 

The following theorem is known as Caratheodory's theorem: 

Theorem 3.3. Given a set of points {xi}™ =l C suppose y G conv{:cj}™ =1 . 
ITien f/iere exists a subset I C {1, ...,n} suc/i t/iat y £ conv{xj}j G 7- and 
{xi}i£i is affinely independent. 

Corollary 3.4. Suppose {ipi}f =i C C d is a scalable frame. Then there is a 
subset {(fi}i € j which is also scalable and is linearly independent. 

Given a unit norm frame {ifi}f =1 C C d we define the set 

n 1 
P({<Pi}i=i) := {(wi,...,w n ) : Wi > 0,^2wi<pi<p* = -I d }. 

i=l 

Proposition 13.21 tells us two things about this set: first we have that w € 
V({^fi}2 = i) if and only if d ■ w is a scaling of {(pi}f =i , and second, that 
V({<Pi}f = i) is a (possibly empty) polytope (see, for example, Theorem 1.1 
in 0). " 

Suppose {ipi}f =i C C d is a scalable frame, and we are given a scaling 
to = w n ). We say the scaling is minimal if : W{ > 0} has no 

proper subset which is scalable. 

Theorem 3.5. Suppose {</?j}" =1 QC d is a scalable, unit norm frame. Ifw = 
(wi, ...,w n ) is a minmal scaling then {(pup* : u>« > 0} is linearly independent. 
Furthermore, V({ipi}™ =1 ) is the convex hull of the minimal scalings, i.e., 
every scaling is a convex combination of minimal scalings. 

Proof. The first statement follows directly from Corollary 13. 41 

We now show that every vertex of V({(pi}f = i) is indeed a minimal scaling. 
Let u E V({(fi}f =1 ) be a vertex and assume to the contrary that u is not 
minimal, then there exists a v £ P such that supp(w) C supp(u). Let 
w(t) = v + t(u — v), and to = mm{ v v J; u : Vi > ui}. We observe that to > 1 
and w(to)i > since supp(v) C supp(u). This means w(to) G P, and u lies 
on the line segment connecting v and w(to) which contradicts the fact that 
u is a vertex. 
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Finally we show that every minimal scaling is a vertex of V ({(pi}^ =1 ) . Sup- 
pose we are given a minimal scaling w which is not a vertex of V({<pi}f = i). 
Then we can write w as a convex combination of vertices, say w = ^tjUi, 
where we know at least two tj's are nonzero, without loss of generality say 
t\ and t%. Since both t\ and t<i are positive and all the entries of v\ and 
V2 are nonnegative, it follows that supp(ui) U supp(t>2) C supp(w), which 
contradicts the fact the w is a minimal scaling. □ 

Theorem 13.51 reduces the problem of understanding the scalings of the 
frame {<£>j}f =1 to that of finding the vertices of the polytope V({ipi}f =1 ). 
Relatvely fast algorithms for doing this are known, see [2]. 



4. When are outer products linearly independent? 

Since most of the results in this paper deal with linear independence of 
the outer products of subsets of our frame vectors we will address this issue 
in this section. It would be nice if there were conditions on a frame {ipi}f =1 
which could guarantee that the set of outer products {<^j(/?*}™ =1 is linearly 
independent, or conversely if knowing that {ifiip*}f =1 is linearly independent 
tells anything about the frame {<pi}™ =1 . One obvious condition is that in 
order for {<Pi<p*}f = i to be linearly independent we must have n < d 2 , and 
when this is satisfied Theorem 12.11 tells us that this will usually be the case. 

Another condition which is easy to prove is that if {v^liLi is linearly 
independent then so is {<fii(p* }™ = i- The converse of this is certainly not true, 
and since we are usually interested in frames for which n > d this condition 
is not very useful. The main idea here is that while the frame vectors live 
in a d-dimensional space the outer products live in a <i 2 -dimensional space, 
so there is much more "room" for them to be linearly independent. 

Given a frame {(/?j}" =1 we define its spark to be the size of its smallest 
linearly dependent subset, more precisely 

spark({(/9j}" =1 ) := min{|7| : {(fi}iei is linearly dependent}. 

Clearly for a frame {ipi}f = i C C d we must have that spark({y?j}" =1 ) < d+ 1, 
if its spark is equal to d + 1 we say it is full spark. For more background on 
full spark frames see [I]. 

Proposition 4.1. Suppose {ip-i}f =1 C C d is a frame with n < 2d — 1. If 
{(fi}f =1 is full spark then {'■fi'-p* }" =1 is linearly independent. 

Proof. Suppose by way of contradiction that {^i}f =1 is full spark but {(/?j</?*}" = 
is linearly dependent. Then we can write an equation of the form 
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with <Zj > for every i £ I, bj > for every j E J, and I fl J = 0. This 
implies that 

span({c^} ie /) = ImQ^a^y*) 

= Im( = sp&n({ipj} je j). 

But since n < 2d — 1 we have either |/| < d — 1 or | J| < d — 1, so this 
contradicts the fact the {^i}" =1 is full spark. □ 

We first remark that the converse of Proposition 14.11 is not true: 

Example 4.2. Let {e±, e 2 , 63} be an orthonormal basis for C 3 and consider the 
frame {ei, e 2 , e%, e\ + e 2 , e 2 + 63}. Clearly this frame is not full spark and yet 
it is easy to verify that {e\e\, e 2 el, e 3 e^, (ei +e 2 )(ei + e 2 )*, (e2 + e3)(e 2 +e 3 )*} 
is linearly independent. 

Next we remark that the assumption n < 2d — 1 is necessary: 

Example 4.3. Let {ei,e 2 } be an orthonormal basis for C 2 and consider the 
frame {ei, e 2 , ei + e 2 , ei — e 2 }. Clearly this frame is full spark but 

eie\ + e 2 e 2 = I 2 = -((ei + e 2 )(ei + e 2 )* + (ei - e 2 )(ei - e 2 )*). 

Finally we remark that with only slight modifications the proof of Pro- 
postion [4.11 can be used to prove the following more general result: 

Proposition 4.4. // spark({</?j}" =1 ) > s then spark({(^j</?|}™ =1 ) >2s — 2. 

Unfortunately, the converse of Proposition 14.41 is still not true. The main 
problem here is that given any three vectors such that no one of them is a 
scalar multiple of another, the corresponding outer products will be linearly 
independent (we leave the proof of this as an exercise). Therefore it is easy 
to make examples (such as Example 14.21 above) of frames that have tiny 
spark, but the corresponding outer products are linearly independent. 

We conclude our discussion of spark by remarking that in [1] it is shown 
that computing the spark of a general frame is NP-hard. Thus, the small 
amount of insight we gain from Proposition 14.41 is of little practical use. 

Another property worth mentioning in this section is known as the comple- 
ment property. A frame {<^j}™ =1 C C d has the complement property if for ev- 
ery / C {1, n} we have either span({</?j}j g /) = C d or span({(/?j}j e /c) = C d . 
We remark that the complement property is usually discussed for frames in 
a real vector space, but for our purposes it is fine to discuss it for frames in 
a complex space. In [3] the complement property was shown to be necessary 
and sufficient to do phaseless reconstruction in the real case. 

If a frame {(fi}™ =l C C d has the complement property then clearly we 
must have n>2d— \ (if not we could partition the frame into two sets each 
of size at most d— 1) and that in this case full spark implies the complement 
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property. If n = 2d — 1 then the complement property is equivalent to full 
spark, but for n > 2d — 1 the complement property is (slightly) weaker. One 
might ask if the complement property tells us anything about the linear 
independence of the outer products, or vice versa. Example 14.21 above is an 
example of a frame which does not have the complement property but the 
outer products are linearly independent, and Example 14.31 is an example of 
a frame that does have the complement property but the outer products are 
linearly dependent. So it seems like the complement property has nothing 
to do with the linear independence of the outer products. 

Given a frame with the complement property we can add any set of vectors 
to it without losing the complement property. Thus it seems natural to ask 
whether every frame with the complement property has a subset of size 
2d—l which is full spark. This also turns out to be not true as the following 
example shows: 

Example 4.5. Consider the frame in Example 14.21 with the vector e± + 
added to it. It is not difficult to verify that this frame does have the com- 
plement property, but no subset of size 5 is full spark. 

We conclude by noting that as in the proof of Proposition 14.11 a set of 
outer products {c/?j<£>*}f =1 is linearly dependent if and only if we have an 
equation of the form 



with en > for every i £ I, bj > for every j £ J, and J n J = 0. This 
is equivalent to {<£>i}?=i having two disjoint subsets, namely {y?j}j e / and 
{<f>j}jeJ, which can be scaled to have the same frame operator. Thus, deter- 
mining whether {<Pi<p*}'i = i is linearly independent is equivalent to solving a 
more difficult scaling problem than the one presented in this paper. 
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